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Abstract 

Nuclei are unique analyzers for the early stage of the space-time development of 
hadronization. DIS at medium energies is especially suitable for this task being sen- 
sitive to hadronization dynamics, since the production length is comparable with the 
nuclear size. This was the driving motivation to propose measurements at HER- 
MES using nuclear targets, and to provide predictions based on a pQCD model of 
hadronization [1]. Now when the first results of the experiment are released [2,3], one 
can compare the predictions with the data. The model successfully describes with no 
adjustment the nuclear effects for various energies, %, pr, and Q 2 , for different flavors 
and different nuclei. It turns out that the main source of nuclear suppression of the 
hadron production rate is attenuation of colorless pre-hadrons in the medium. An al- 
ternative model [4] is based upon an ad hoc assumption that the colorless pre-hadron 
is produced outside the nucleus. This model has apparent problems attempting to 
explain certain features of the results from HERMES. A good understanding of the 
hadronization dynamics is important for proper interpretation of the strong suppres- 
sion of high-p-r hadrons observed in heavy ion collisions at RHIC. We demonstrate 
that the production length is even shorter in this case and keeps contracting with 
rising pr- 
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1 Introduction 



Recent measurements by the HERMES collaboration [2] of semi-inclusive production 
of hadrons in deep- inelastic scattering (DIS) off nuclei have provided precious informa- 
tion about the space-time development of hadronization. These measurements were 
proposed [1] back in 1995 and predictions were made within a model of hadronization 
based on perturbative QCD. One of the goals was to reach a better understanding of 
in-medium hadronization in order to provide a more reliable interpretation of high-pr 
hadron production in heavy ion collisions considered as a probe for the dense matter 
produced. 

Now, when some of the results of these measurements are released, it is proper to 
compare the predictions with the data and draw conclusions. In the present paper we 
revisit the problem of in-medium hadronization, re-introduce the model, and perform 
calculations for HERMES kinematics. We demonstrate that absorption of the pro- 
duced colorless pre-hadron is the main source of nuclear suppression. Our predictions 
need no adjustment and are in good agreement with the data. 

Older models [5-9] were based on the string model which were quite pedagogical 
and helpful for intuitive understanding of the space-time development of hadronization. 
On the other hand, one could not consider such a nonperturbative phenomenology as 
a realistic scheme for hard reactions like DIS. These models had a very little predictive 
power, no access to Q 2 and pr-dependences, missed the effect of color transparency, 
etc. 

The release of data from the HERMES experiment has stimulated a new wave of 
theoretical models. All of them have low predictive power and are fitted to data which 
are supposed to be explained. This is why all of them agree with data, although they 
explore quite different physical ideas. 

Some of the models [10-12] continue developing already known ideas based on the 
phenomenology of nonperturbative hadronization. Others [4,13], employ the idea of 
perturbative induced energy loss, but push to the extreme. Namely, they make two 
strong ad hoc assumptions which have no justification. First, it is assumed that pre- 
hadrons are always produced outside the nucleus. Second, the effect of induced energy 
loss is accounted for via a simple shift of the variable in the fragmentation function, 
which is equivalent to the assumption that hadronization starts only after the leading 
parton leaves the nucleus. 

Since the induced energy loss scenario is a popular model or the nuclear suppression 
of high-pT events observed at RHIC, we confront some of the predictions of this model 
with data from HERMES. We sort out those observables which are sensitive to the 
model assumptions. 

1.1 Why nuclear target? 

Particles produced in hadronic collisions and detected at macroscopic distances carry 
limited information about the hadronization dynamics. The most important details 
are hidden at the early stage. Nuclear targets provide a unique opportunity to look 
at the early stage of hadronization at distances of a few Fermi from the origin. A 
quark-gluon system originated from DIS propagates through the nuclear medium and 
interacts with other nucleons. Modification of the differential cross section of particle 
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production can bring forth precious information about the structure of the excited 
system and its space-time development. 

We assume in what follows that Bjorken x is sufficiently large, x ~ 0.1, therefore the 
lepton interacts incoherently with only one bound nucleon, no shadowing is possible. 
Besides, this region is dominated by valence quarks, i.e. one can treat the DIS as 
electron-quark scattering with energy transfer, v, to the knocked out quark. At small 
x dominated by the sea the process of DIS looks different (in the nuclear rest frame) : 
the virtual photon produces two jets, q and q, which share the full energy v. Therefore, 
the variable Zh, which is the fraction of v transferred to the detected hadron, cannot 
reach the kinematic limit z^ = 1 for either of the two jets, unless one of them takes 
the whole energy of the photon. This fact makes the ratio of nuclear to nucleon cross 
sections [see Eq. (5)] fall faster at Zh — ► 1 than at large x. Additionally, negative kaons 
which at large x are produced mainly from gluons, at small x are generated by the 
same mechanism as positive kaons. Thus, one should be cautious about the range of 
x involved in the analysis, and make a proper x-binning of DIS data. 

One can also use the process of hadronization as a tool for probing the medium 
properties. This is the driving idea of the so called "jet quenching" probe for creation 
of dense matter in relativistic heavy ion collisions [14]. This tool, however, can work 
properly only if the hadronization dynamics is reliably understood, which seems to be 
still a challenge. 

1.2 Vacuum and induced energy losses 

Energy loss is a hot topic nowadays, however, related discussions are confusing some- 
times, since vacuum and induced energy losses are mixed up. Due to confinement, an 
energetic parton cannot propagate in vacuum as a free particle, but hadronizes and 
produces a jet of hadrons. The parton shares its energy with the produced hadrons, 
therefore it is gradually losing energy. We call this vacuum energy loss. As a guidance 
from string model [15] the rate of energy loss should be constant, of the order of the 
string tension n = lGeV/fm (see Sec. 2). Therefore, the energy loss rises linearly 
with time or length of the path, 

AE vac (L) = kL . (1) 

If the parton originates from a hard reaction (DIS, high-p^ process, etc), gluon 
bremsstrahlung should be an additional or even the main source of energy loss. Indeed, 
as a result of a strong kick the parton should shake off a part of its color field with 
transverse frequencies controlled by the strength of the kick. Apparently, energy loss 
may be very extensive and may substantially exceed the static value given by the 
string model. In this case an approach based on perturbative QCD should be more 
appropriate. 

One can trace similarity with energy loss in QED. If a plate of a big capacitor, 
which transverse dimensions are much larger than the distance between the plates, is 
accelerated, the main contribution to the retarding force comes from the static field. 
Indeed, the radiation from the edges of the capacitor is small and can be neglected. 
This an analog to a color string (tube) model which might be a reasonable approxima- 
tion for soft processes. Indeed, this case is associated with constituent quarks which 
size is of the order of the tube length. 
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On the contrary, if one accelerates point-like charges, the main contribution to 
energy loss should be bremsstrahlung, rather than Coulomb forces. This would be an 
analogy for quark jets originating from DIS. 

Amazingly, the radiative energy loss turns out to rise linearly with the path length 
[16] like in the string model, 

AE vac (L) = ^-a s (Q 2 )Q 2 L . (2) 

This happens because the gluons lose coherence with the source, i.e. are radiated at 
different time intervals. 

In the case of hadronization inside a medium, an additional source of energy loss 
is interaction of partons with the medium. Due to multiple collisions the parton 
increases its transverse momentum squared linearly with the path length, since it 
performs Brownian motion in the transverse momentum plane. Experimentally, the 
broadening of transverse momentum is a rather small effect. For instance, the broad- 
ening measured in the Drell-Yan reaction on a nucleus as heavy as tungsten is only 
Ap\ « 0.1 GeV 2 . Apparently, the induced energy loss [18] generated by such a small 
momentum transfer is a small addition to the vacuum energy loss originated from a 
hard process (at least in cold nuclear matter), 

AE ind {L) = ^a s Ap 2 T L = ^ a s C(E) p A L 2 , (3) 
where C{E) = daq q (rT , E) / dr\ [19,20] is the derivative of the universal phe- 

Tx =0 

nomenological dipole cross section which depends on the transverse dipole separation 
tt] pa is the nuclear density. Of course, application of perturbative QCD in this case 
is questionable. 

The intuitive interpretation of these results is rather straightforward. It takes a 
time, called the coherence time, to radiate a photon or gluon which becomes incoherent 
with the source [1], 

(4) 

f\j>j-i 

where kx and a are the transverse momentum and fraction of energy E taken away 
by the radiated quantum. According to the Landau-Pomeranchuk principle, radiation 
at long coherence times does not resolve the details of interaction, whether it was 
a single kick, or a few of them. What matters is the final transverse momentum of 
the parton. Since the additional transverse momentum gained by the parton due to 
multiple interactions in the nuclear medium is very small, the corresponding correction 
(induced) to the energy loss is small as well. 

This is why it is not easy to see in data the effects of induced energy loss. In 
particular, the analysis of nuclear effects in Drell-Yan reaction performed in [21] pro- 
vided information about vacuum energy loss. Indeed, according to the kinematics of 
the E772/E866 experiments (and other experiments at lower energies) most of events 
correspond to a short coherence time for the Drell-Yan process which occurs nearly 
momentarily deep inside the nucleus. Prior to that the incoming hadron interacts 
softly on the front surface of the nucleus. This collision breaks down the coherence 
of the projectile partons which start losing energy on the path from the soft collision 
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point until the point where the Drell-Yan reaction takes place. One should expect the 
vacuum value for the rate of energy loss which is at least —dE/dz « lGeV/ fm or 
higher, much larger than what one could expect for induced gluon bremsstrahlung. In- 
deed, the analyses [21] led to —dE/dz = (2.8 ±0.4 ±0.5) GeV/ fm, which is a summed 
effect of vacuum and induced energy losses. 

Naively, one might expect that the effects of vacuum energy loss cancel out in the 
nucleus to nucleon ratio, since they are identical in the numerator and denominator. 
This is not correct, firstly, induced energy loss is a part of hadronization and it does 
not result in a simple shift of the variable in the fragmentation function. Secondly, 
the vacuum energy loss controls the time scale of hadron production which has a great 
impact on nuclear effects (see next section). 

1.3 Production and formation times 

One should discriminate between formation of the final hadron and production of a 
colorless state which is not yet an eigenstate of the mass matrix and may be projected 
to various hadronic wave functions. We call such states pre-hadrons and the related 
time scale production time, t p . The properties of such pre-hadron states and their 
attenuation in nuclei have been intensively studied over the last two decades, both 
theoretically and experimentally (see e.g. in [22]). It has been proven that such 
pre-hadrons of a reduced size are indeed produced, and that nuclei are indeed more 
transparent for such states. A signal of color transparency in diffractive DIS was 
observed recently in the HERMES experiment [23,24]. 

To form the hadronic wave function the constituents of the pre-hadron have to 
circle at least once along their orbits. Dilated by Lorentz transformation this time 
scale, which we call formation time, tf, is proportional to Zh- 

On the other hand, the production time t p was proven in [6, 25] and confirmed 
within the LUND model in [26] to vanish proportionally to 1 — Zh at Zh — ► 1. The 
latter property follows from energy conservation: the longer the hadronization process 
is lasting, the more energy is lost by the leading quark, unless its energy falls below 
the desired hadron energy. In particular, at Zh = 1 no energy loss for hadronization is 
allowed, so the colorless ejectile has to be created instantaneously. 

Note that an averaged hadronization process is lasting a long time proportional 
to the energy of the parton initiated the jet. However, the mean energy fractions z^ 
of the produced partons are quite small. Selecting rare events with a leading hadron 
produced with large Zh one imposes a severe restriction upon the hadron production 
time which must be proportional to (1 — z^j- 

1.4 What are the observables? 

The nuclear medium affects the momentum distribution of the produced particles, 
which is called sometimes a modified fragmentation function Df^{zh,pr-, Q 2 ,v), where 
Pt is the transverse momentum of the hadron, and v and Q 2 are the energy and 
virtuality of the photon respectively. The energy dependence of signals that 

QCD factorization is broken and D^ff cannot be treated as a fragmentation function, 
e.g. the QCD evolution equations cannot be applied. Only at very large v and Q 2 
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nuclear effects disappear and factorization is restored. In this trivial limit, however, 
all physical information we are interested in is missed. 

Experimental results are usually presented in a simple form, as a ratio of the 
nucleus-to-nucleon hadron multiplicities, 

2 dnjyA -> hX)/dz h d 2 p T 

R A (z h , PT ,Q ,„) = Adn{rN ^ hx)/dzhd2pT ■ (5) 

Because of limited statistics this four-dimensional ratio is usually presented as one- 
dimensional ratios integrated over other variables. Some of correlations are, however, 
very informative. For instance, the pt dependence of the ratio binned in brings forth 
precious information about the space-time pattern of hadronization (see Sec. 9.3). 

Nuclear modification of the hadronic spectra was considered for high-pj- hadron 
production in [6,27], for DIS in [5,7,9] and for hadroproduction of leading particles 
on nuclei in [26, 28]. 



1.5 What can we learn? 

The questions, which can be answered in such an analysis are: 

- How long does it take to produce a pre-hadron? 

Nuclear suppression of hadron production, especially its Zh and pr dependences, are 
sensitive to the production time t p . 

- How does the produced pre-hadron evolve and attenuate in nuclear matter? 

Since we are interested in production of leading particles, we should forbid inelastic 
interaction of the produced colorless wave packet at times t > t p . It would lead to a new 
hadronization process which ends up with a smaller value of Zh- The corresponding 
correction is easy to estimate and it turns out to be small. The condition of no 
interaction results in an attenuation of the produced wave packet in nuclear matter, but 
the absorption cross section of the pre-hadron may be different from the hadronic one. 
Particularly, in DIS with high photon virtuality Q 2 the produced wave packet may have 
a smaller transverse size, similar to what is known for exclusive particle production. 
Correspondingly, the color transparency effects [29, 30] should be important and we 
take care of that in what follows. 

- Does a fast quark attenuate in nuclear matter? 

At first glance an energetic quark should not attenuate since it cannot be stopped 
or absorbed in the medium. Multiple soft interactions can only rotate the quark in 
color space, while the energy loss is a small fraction of the quark energy. However, by 
"attenuation" of a quark, we mean the suppression of the production rate of the final 
hadron produced with a certain momentum, related to the quark interactions in the 
nuclear medium. In the string model any color-exchange interaction leads to formation 
of a new string, but with reduced initial energy. This fact results in a reduction of 
hadron production at large Zh — > 1 [8]. Moreover, the leading quark may pick up new 
strings due to presence of higher Fock components, giving rise to induced energy loss 
in this model [21]. If one tries to treat the reduced production rate of hadrons as an 
effective attenuation of the quark in the medium, the effective quark absorption cross 
section cannot be treated as an universal parameter. Its value essentially varies with 
geometry and kinematics [8]. In the present paper we include induced energy loss 
explicitly. 
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2 Lessons from the string model 



Apparently, the nonperturbative string model cannot be a realistic approach to such 
a hard process as DIS. The model does not have any access to (5 2 -dependence, color 
transparency, etc., which are just the phenomena that are crucial for understanding the 
properties of in-medium hadronization. At the same time, the model is simple, very 
intuitive, and helps to understand the general features of the space-time development 
of hadronization. 

As was already mentioned, we assumed that data are taken at large Bjorken x > 0.1 
to make sure that valence quarks in the target nucleon give the main contribution. A 
color string 1 formed between the quark and the debris of the target nucleon, slowing 
down the former and speeding up the latter. Naively, one may think about a long string 
stretched across the nucleus. Instead, a very short object is propagating through the 
nucleus. Indeed, a simple kinematics shows that the maximal possible length of the 
string in the rest frame of the nucleus is, 

Lmax = i (6) 
K 

where m is the mass of the debris of the target nucleon (e.g. a diquark) attached to the 
slow end of the string. The string tension k = l/(2ira' R ) ~ 1 GeV/fm is related to the 
slope a' R of the Regge trajectories [15]. Since the nucleon debris are probably lighter 
than the nucleon, L max < 1 fm. In reality the string should be even much shorter due 
to the Schwinger phenomenon, spontaneous qq pair production from vacuum breaking 
the string to shorter pieces. 



Thus, a heavy but short string with the effective mass M ~ yfs = y m? N — Q 2 + 2m^u 
propagates through the nucleus. Due to qq pair production in the color field of the 
string light pieces of the string are chipped off, while the leading heavy part of the 
string keeps propagating on. Apparently, the leading quark is slowing down, losing 
energy with a constant rate, 

dE 

where E(t) is the quark energy and t = z are the time and longitudinal coordinate. 
Note that the rate of energy loss (7) is a constant which is independent of the quark 
flavor, energy and virtuality, and is invariant relative to longitudinal Lorentz boosts. 
Correspondingly, the total energy loss rises linearly with the length Az of the path, 
AE = kAz. The lost energy goes to acceleration of the target debris and production 
of new hadrons. The process of hadronization completes when the mass of the leading 
piece of the string is reduced by many decays down to the hadronic mass scale and 
the leading pre-hadron (or a cluster) is produced. 

How does the production time depend on the initial quark energy v and its fraction 
Zh carried by the produced hadron? For low energy hadrons it is short, but rises with 
their energy, t p f» vz^jlK. One might think that such a linear 2^-dependence is an 
unquestionable truth since it is dictated by Lorentz time dilation. It was claimed, 
however, in [6,25] and confirmed within Lund model [26], that there is an opposite 
trend for leading hadrons, namely, the production time vanishes like t p = (1 — Zh)v / k 



: Of course a one-dimensional string is an idealization. We neglect the transverse size of the color- flux 
tube [15], unless otherwise specified. 
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at Zh — > 1. Although it was first derived in the string model [6,25], this is a general 
property, since is dictated by conservation of energy. Indeed, unless color neutraliza- 
tion happens, i.e. a colorless system with energy ZhV is produced, the leading quark 
keeps losing energy. Therefore, the smaller (1 — z^) is, the faster this process must 
be completed, otherwise the energy needed for the leading hadron production will be 
wasted. 

Note, that the produced colorless object should be treated as a pre-meson. It may 
take a long time proportional to z^u to form the final hadron wave function. We call 
this formation time, and it is always longer than the production time, tf > t p . 



Thus, we arrive at two different types of end-point behavior of the production 
time [6,25,26], 



The concept of production time is crucial in the case of a nuclear target: a colorless 
pre-hadron produced inside the nucleus should not interact (inelastically) on the way 
out of the nucleus, otherwise hadronization will be triggered, and the final hadron will 
emerge with a substantially reduced energy. This restriction leads to an attenuation 
which depends on the interaction cross section of the pre-hadron (see below) . For this 
reason different processes on nuclei are suppressed near the kinematic limit z^—^l, for 
example, back-to-back high-pr dihadron production, DIS, forward particle production 
in hadron-nucleus interactions, etc. 

As was mentioned above, the production time Eq. (8) corresponds to creation 
of a fast colorless piece of the string, a pre-hadron, rather than the final hadron. 
The latter is supposed to have a wave function describing the specific distribution of 
quark momenta inside the hadron, different from the quark momentum ordering in 
the string. The wave function of the hadron is developing over a longer period of 
time (formation time), which is proportional to the hadron energy, z^is, rather than to 
(1 — Zh)v. These two time scales are mixed up sometimes. At low Q 2 the pre-hadron 
has a large absorption cross section, and nuclear attenuation is controlled mainly by 
the production time Eq. (8). In this case the formation time has a moderate influence 
on nuclear effects. On the contrary, at large Q 2 the formation time scale is more 
important for nuclear attenuation. We will include both effects in our calculations. 

Another question which may be answered within the string model is how multi- 
ple interactions of the hadronizing quark at earlier stage, t < t p affects the hadron 
production rate? Sometimes [2, 7] this attenuation is treated as absorption with an 
effective quark-nucleon cross section treated as a universal fitted parameter. However, 
the physical origin of such an attenuation discussed in [8] shows that it cannot be de- 
scribed in terms of a universal parameter like a quark absorption cross section. Indeed, 
every interaction of the leading quark leads to a new hadronization process starting 
from the very beginning, but with a reduced initial energy E q (z) = v — kz, where z is 
the distance covered by the quark. Since the energy of the final hadron is fixed, the 
variable of the fragmentation function should be redefined as, z^ => z^/(l — kz/u). 
Apparently, the increase of Zh should lead to a suppression which depends on v and 
Zh and cannot be treated as a universal parameter [8]. 

Moreover, multiple quark interactions at t < t p induce an additional energy loss 




z h < 1 
1 - z h < 1 



(8) 
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with a rate rising linearly with the path length z [21], 

(9) 

where (n(z)) = p^z is the mean number of collisions experienced by the quark 
over the distance z; af^ ls the model dependent effective quark-nucleon inelastic cross 
section. The first term in (9), same as in (7), is the usual retarding force produced 
by the string, while the second term is due to additional strings created by multiple 
interaction of the quark. These extra terms correspond to higher Fock component of 
the quark and can be motivated within the dual parton model [31] (more details can 
be found in [21]). 

Thus, the total energy loss on a distance L contains a correction oc I? . 

AE(L) = kL[1 + \af lPA L] . (10) 

In the additive quark model of^ ~ 10 mb, which is the same value as follows from the 
dual parton model [31]. Note, that this value is twice as large as in the more realistic 
model [32] , which demonstrates that the dominant part of the hadronic cross section 
does not obey quark additivity. 

It is questionable whether one can apply this phenomenology to DIS. At this point 
we feel that our unjustified assumption that soft interaction phenomenology, like the 
string model, is relevant for final state interaction in DIS has been pushed to the 
extreme, and we should consider a more realistic approach based on perturbative 
QCD. 



a Hi 

H 7. 



= —K 



1 + (n(z)) 



3 Radiative energy loss: hadronization in 
vacuum 

3.1 Vacuum energy loss 

When a high-momentum parton is produced in a hard reaction, it triggers gluon 
bremsstrahlung, which lasts until the moment of color neutralization. Gluons are not 
radiated instantaneously, but it takes time to lose coherence with the parent parton. 
This quantum-mechanical uncertainty called coherence time can be introduced differ- 
ently. A strict definition comes within precise calculations and is related to interference 
between quanta radiated at different positions. For instance, two amplitudes of gluon 
radiation by a quark at different coordinates separated by the longitudinal distance 
Az get phase shift Azqi, where qi = (Mq — m q i )/2E q is a difference between the 
longitudinal momenta of the quark and the produced quark-gluon final state which 
has effective mass squared, 



k\ + am q 
n,: = a(l-a) ' 



Mk = ' ■ (11) 



Here kj- and a are the transverse momentum and fraction of the light-cone momentum 
of the quark carried by the gluon; E q = v is the initial quark energy (we assume that 
Bjorken x is large and the quark gets the whole energy transfer). 
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Apparently, only those gluons can interfere, which are radiated within the coher- 
ence time interval, or coherence length (we assume z = t and neglect the quark mass 
hereafter), given by Eq. (4) with E = v. This time may be also interpreted as the 
lifetime of the quark-gluon fluctuation, or as the formation time of radiation. At any 
rate, only those gluons may be treated as radiated during time t, which have coherence 
time t c < t, otherwise they are still coherent with the quark and should be considered 
as a part of its color field. 

Correspondingly, one can estimate the energy deposited into the radiation as 

Q' 2 1 

AE(t)=u f dk\ ( daa-^-&{t-t c ) (12) 

A2 

where A is an infrared cutoff which we fix at A = Aqcd = 200 MeV. The number of 
radiated gluons dnc/dadk^ = j/ak 2 ^, where 7 = 4a s (k^) /3n [16,33]. As usual, we 
employ the approximation of independent and soft radiation, a <C 1. 
Eq. (12) leads to a linear time dependence of energy loss, 

AE{t) = 1 -{Q 2 -\ 2 )t (13) 

i.e. the rate of energy loss is constant, as it is for strings. Although gluons are 
continuously radiated with different coherent times and different energies, the amount 
of energy taken away per unit of time remains the same. This interesting observation 
was demonstrated in [16,17]. 

The radiated energy is distributed among the produced hadrons, but we are inter- 
ested in the most energetic one, the leading hadron, which carries the largest fraction 
Zh of the initial momentum of the quark. This condition and energy conservation 
impose additional restrictions on the integration in Eq. (12), resulting in a different 
time development of the hadronization and the energy loss. 

If the leading hadron is produced with large Zh and includes the leading quark, 
none of the radiated gluon may leave with a fraction of the initial quark momentum 
exceeding a > 1 — z^. For instance, in the limiting case Zf- L — > 1 any radiation is 
forbidden, which leads to a so-called Sudakov's suppression of the cross section. With 
decreasing Zh, this restriction gradually softens, but still should not be forgotten. 

In order to estimate the time-dependence of radiative energy loss specified for the 
case of the leading hadron production we introduce energy conservation in Eq. (12) 
via a step-function 0(1 — Zh — a). The result of integration reads, 

AE(t) = It (Q 2 - A 2 ) 9(ti - t) 

+ {1/7(1 1 + ln {j;) ~ \ Xh } (* - *i) e (*2 - <) 

+ n (l- Zh )ln(^pjQ(t-t 2 ) (14) 



where 



2u , 1- z h 

Q z xitin 
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*2 = ^2*1 , (16) 

and x = Q 2 /2mNU is Bjorken variable. 

The time dependence of radiative energy loss is illustrated in Fig. 1 assuming 
Q 2 3> A 2 . During the time interval t < t±, the rate of energy loss is constant, —dE/dt = 
7Q 2 /2, exactly as in the case of Eqs. (12), (13) with no restriction on the radiated 
energy. 



AE(t) 




ti t 2 
t 

Figure 1: Time-dependence of energy loss corresponding to Eq. (14). The 
notations are explained in the text. 



Over a longer time interval, more energetic gluons can be emitted in accordance 
with Eq. (4) and the restriction a < 1 — zy t becomes effective. As a result, the loss of 
energy slows down to a logarithmic t-dependence. This implies that the rate of energy 
loss decreases as function of time as dE/dt = —7^(1 — Zh)/t. 

3.2 Sudakov's suppression 

According to Eq. (14) at much longer times either very energetic gluons can be emitted 
(which may be forbidden by the condition a < 1 — z^), or gluons with very small kx- 
The latter are cut off in Eq. (12). This completely eliminates any energy loss at long 
times t > t2 as one can see from Eq. (14) (see also Fig. 1). This sounds puzzling, unless 
the suppression due to a Sudakov-type factor is taken into account. This suppression 
becomes active t > t±, when the production of the leading hadron with large Zh 
restricts the spectrum of radiation. 

Summing up all the nonradiated gluons we get the following Sudakov's suppression 
factor, 

S(t,z h ,Q 2 ,v) = exp -h G (t, z h ,Q 2 ,v) , (17) 

which corresponds to the Poisson distribution of independently radiated soft gluons. 
The number of gluons which have sufficiently short coherence time to be radiated 
during time interval i, but conflicting with the leading hadron production since they 
are too energetic, a > 1 — Zh, is given by, 

1 Q 2 2 

n G (t,z h ,Q 2 ,v)= 7 J ^ j ^jqi-exp(-t/t c )] . (18) 
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To be more realistic we smoothed out the step- function 0(t — t c ) here, as well as in 
(14). We assumed that the radiation rate of gluons with coherence time t c has a decay 
rate distribution, 

= !.-«*■. (i9) 

dt t c v ' 

Then the probability to radiate a gluon during time interval t is, 

P(t) = 1 - e~ t/tc , (20) 

which is a replacement for the step- function 0(i — t c ) in (14) and (18). 

The Sudakov formfactor, Eq. (17), substantially reduces the mean time interval of 
hadronization and the production time of the leading pre-hadron. 



3.3 Leading hadron production 

The estimate Eq. (1) for energy loss was obtained within the string model, which 
has no access to Q 2 dependence and assumes that the final hadron is immediately 
produced once the quark energy degrades down to z^v. Perturbative QCD gives 
the opportunity to develop a more detailed and profound model for hadronization 
and color neutralization. In the large iV c limit each radiated gluon is equivalent to 
a qq pair, and the gluon bremsstrahlung can be seen as production of a system of 
colorless dipoles (compare with [34,35]). This would be a nonperturbative mechanism 
of hadronization. It is more appropriate to treat hadronization, i.e. the transition 
G — > qq, perturbatively if the production time is short. Then it is natural to assume 
that the leading hadron originates from the color dipole which includes the leading 
quark and the antiquark from the last radiated gluon. In the large N c approximation 
they are automatically tuned in color, i.e. this pair is colorless. This qq dipole is to 
be projected on the hadron wave function, ^fh(P, It), where (3 is the fractions of the 
hadron light-cone momentum carried by one of the quarks, and It is the transverse 
momentum of the quarks in the hadron. 

Let us evaluate the distribution function W(t, Zh,Q 2 , v) for production time t of a 
pre-meson, 



W(t,z h ,Q 2 ,u) = N [ —5 
J a 



Zh 



Q 

I 

A 2 
1 



dk T exp(— t/t c 



x J dp 5 
o 



P- 



a 



2-a 



I 



dl 2 5 



ij, — j^k T 



\* h (P,l T )\ 2 S{z h ,t,Q 2 ,u) . 



(21) 



Here a is the fraction of the quark light-cone momentum carried by the gluon emitted 
at the time t by the quark of energy E q (t) = v — AE(t). The mean radiation time 
of this gluon t c is given by Eq. (4). Assuming that the radiated gluon splits to a qq 
pair which shares the gluon momentum in equal parts, we conclude that the leading 
colorless qq pair is produced with an energy E q (t)(l — a/2). On the other hand, this 
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must be the observed energy z^u of the produced hadron, the condition controlled by 
the first 5-functions in Eq. (21). The transverse momentum It of the quark and the 
antiquark within the colorless dipole is related to for as It = 3kT/^, while the fraction 
of the hadron light-cone momentum carried by the antiquark is f3{t) = a(t)/[2 — a(t)]. 
This is why we need the last two (^-functions in Eq. (21). The normalization factor ./V 
in (21) cancels in ratio (5), so we do not specify it. 

Note that Eq. (21) has a probabilistic structure; instead of projecting the pro- 
duction amplitude on the hadron wave function, we convolute their squares. This 
simplification is related to our need to include the time development of hadroniza- 
tion and make use of the coherence time. To do it correctly one should work with 
amplitudes taking care of interferences at all stages. In principle, this can be done 
within the path integral approach (partially applied below), however, technically it 
is so complicated that the exact solution of this problem does not look feasible at 
present. Nevertheless, we will see in what follows that the approximation Eq. (21) 
does a good job describing data. 

We use the parameterization for the hadronic wave function in (21) in the form of 
the asymptotic light-cone meson wave function [24], 



OC 



0(1-/3) 
0(1 -(3)+a 



exp 



0(1-/3) + ao 



(22) 



where R\ = 8(r 2 ) em /3 is related to the mean hadronic electromagnetic radius squared, 
{ r h)em- Normalization is not important since it can be absorbed into the factor N 
in (21). The parameter ao in (22) incorporates confinement for highly asymmetric 
qq configurations with — > 0, 1. We fix this parameter demanding the mean charge 
radius of the hadron to be correctly reproduced, 



(P 2 ) 



r d§ 

J /3 (1-13) 



fdp 2 \* h (f3,p)\ 2 p 2 8 



= -(r z ) 

g V / em ' 



(23) 



where ^(0, /?) is the hadronic wave function in coordinate representation, i.e. the 
Fourier transformation of Eq. (22). For production of pions with (r 2 ) em = 0.44/m 2 [36] 
we found ao = 1/12. 

Examples for the production-time distribution W(t) at < v >= 12GeV, < Q 2 >= 
3GeV 2 and different are presented in Fig. 2 2 . We see that as a function of z^ 
the mean production time decreases similar to the expectations of the string model 
(8) (but with different value of k), which is not a surprise, since this behavior reflects 
conservation of energy. 

There are no data, of course, for W(t, Zh, Q 2 ) to compare with, but the integral over 
t at fixed Zh is the probability to produce the hadron with energy z^v. This probability 
is the fragmentation function which is rather well known phenomenologically: 



D h/q (z h ,Q 2 ) = Jdt W(t,z h ,Q 2 ) . 



(24) 



Tilde here means that this function describes fragmentation of a quark only to a 
(leading) hadron which includes the parent quark. 



2 Note that a function, t W(t) was plotted in Fig. 1 of Ref. [1], but a factor t was missed in the notations. 
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Now we can fix the normalization factor N in (21). In fact, the normalization 
of Dfr/g is different from that for the conventional fragmentation function which is 
normalized to the particle multiplicity, Jq 1 dzh D h / q (zh,Q 2 ) = (n/J. Apparently, only 
one hadron in the produced jet contains the original quark, therefore, 

1 

J dz h D h/q (z h ,Q 2 ) = l. (25) 
o 

Now we are in a position to compare the modeled fragmentation function with a 
realistic phenomenological one. In view of the specifics of our model which is designed 
to reproduce production of leading hadrons, we may expect agreement only at large Zh, 
probably z^ > 0.5. Our calculations for Q 2 dependence of D h / q (zh, Q 2 ) at different z^ 
are depicted in Fig. 3, and are compared with popular parametrizations fitted to data, 
[37] and [38] shown by thin and thick bars respectively. The difference between the 
parametrizations may be treated as a systematic uncertainty. The model reproduces 
the phenomenological fragmentation function rather well. This encourages us to apply 
this approach to nuclear targets as well. 

4 Evolution and attenuation of the pre-hadron 
in the nuclear environment 

Production of a leading colorless pre-hadron with the detected momentum completes 
the process of hadronization. Further propagation of the pre-hadron in vacuum and 
formation of the hadronic wave function have no observable effect, since it leads only to 
a phase shift. This changes drastically in the presence of an absorptive medium. In this 
case the production time t p is of special importance. Indeed, when color neutralization 
happens and a pre-hadron with the desired energy, Zhf, is produced, any subsequent 
inelastic interaction in the medium gives rise to a new hadronization process, and the 
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Figure 3: The model prediction for the fragmentation function D(zh, Q 2 ) 
in comparison with parametrizations fitted to data, [37] (thin vertical 
bars) and [38] (thick bars). 



quark energy substantially degrades down to a lower value of Zh- This means that 
the pre-hadron should have been produced with a larger z' h > z^, i.e. with a smaller 
probability. The suppression of the hadron production rate turns out to be so strong 
that we can neglect these double-step processes. 

If the energy of the pre-hadron (qq dipole) is sufficiently high to freeze the trans- 
verse qq separation (due to Lorentz time dilation) for the time of propagation through 
the medium, then the survival probability of the pre-hadron, which we also call nuclear 
transparency, gets the simple eikonal form [39], 



Tr 



^h{r T ) 



exp 



2 u qq\ 



(r T )T A ]\^g q (r T )) 



<*ft(rr)|tf«(rr)} 



(26) 



where Ta is already introduced nuclear thickness, given in this case by integral of the 
medium density along the path of the dipole, Ta(z\,Z2) = dz p A {z). The universal 
dipole-nucleon cross section, cr^(r-r), was introduced earlier, '^(r-r) and ^(rr) are 
the light-cone wave functions of the initial qq-dipole (pre-hadron) and the final hadron 
respectively. 

In the case of medium energies (HERMES, Jefferson Lab) the dipole size may sub- 
stantially fluctuate during propagation through the nucleus and nuclear transparency 
might substantially differ from (26). The rigorous quantum- mechanical approach in- 
corporating all effects of fluctuations should include all possible trajectories of the 
quark and antiquark, rather than just fixed quasi-classical ones. This path integral 
technique leads to an expression for the survival probability of a dipole in a medium 
similar to (26). However, the exponential attenuation factor in the numerator must 
be replaced by the light-cone Green function, Gq q (z2- I ri \ z±,fi), which describes prop- 
agation of the qq pair with initial and final separations r\ and r*2, between points with 
longitudinal coordinates z\ and Z2 [24,39,40], as is illustrated in Fig. 4. Then the 
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Figure 4: The light-cone Green function Gq q (z 2 , f 2 ', zi, fi) describing 
propagation of a qq pair over all possible paths in a medium between 
points with coordinates Z\ and z 2 . 



nuclear transparency Eq. (26) takes the form, 



Tr( Zl ,z 2 ) = 



J d 2 nd 2 r 2 y* h (r 2 )G qq 



[z 2 ,f 2 ;zi,ri) i & 



qq 



Jd 2 r^* h (r)^q q (r) 



(27) 



The Green function satisfies the two-dimensional light-cone Schrodinger equation, 



G ijq ( Z2 ,f 2 ;z 1 ,f 1 ) 

dz 2 



A 



2E h (3(l-(3) 



+ V qq (z 2 ,r 2 ,P) 



x Gg q (z 2 ,f 2 ; zi,n) , 
6W( 



(28) 



with the boundary condition Gq q (z 2 , r 2 ; zi,r\)\ z , 2=zi = 5^ ' (ri — r 2 ). Here Eh = Zh v is 
the hadron energy, and (3 is the fraction of the hadron light-cone momentum carried 
by the quark. The first, kinetic, term in squared brackets is responsible for the phase 
shifts which are related to quark masses, m q , and their transverse motion described 
by Laplacian A r acting on coordinate r. 

The imaginary part of the light-cone potential in (28), 



lmVqq(z,f,p) = --(Jq q {r) p A {z) , 



(29) 



is responsible for attenuation of the qq in the medium, while the real part represents 
the interaction between the q and q. The latter is supposed to provide the correct 
light-cone wave function of the final meson. For the sake of simplicity we use the 
oscillator form of the potential, 



ReVgq(z,f,f3) 



a 4 ((3)r- 



2 E h 0(1-/3) 



(30) 



which leads to a Gaussian r-dependence of the light-cone wave function of the meson 
ground state. The shape of the function a((3) was discussed in [41], and fixed by 
photoproduction data a 2 ((3) = Oq + a 2 (3(1 - /?), where Oq = v 1 - 15 x (112 MeV) 2 ; a 2 = 
(1 — -u) 1 - 15 x (165 MeV) 2 , and v takes any value in the interval < v < 1. 

If we assume also that the small-r dependence for the dipole cross section, cr^(r) = 
C(Eh) r 2 , then Eq. (28) has an analytic solution [42], 



G, 



qq{Z2,r 2 ;z 1 ,r 1/ 



a 2 (f3) 



2 7r i sin(^ Az) 
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x exp 
x exp 

where Az = Z2 — Z\ and 



\ 2sin(£ Az 



i mi Az 



{r\ + r|) cos(£ Az) — 2 r \ ■ r 2 
2E h p(l-~Jj) 

yV(/3) - jgg - P)E h C(E h ) p~ A 
#h 0(1 - 0) 



(31) 



(32) 



This solution is valid only for a medium with a constant density. However, using 
the recurrent relations between Green functions for different space intervals one can 
analytically solve the problem for a varying density as well [39]. 



5 Nuclear attenuation caused by absorption 

We assume that the colorless qq pair (pre-hadron) produced at distance l p from the DIS 
point has a Gaussian distribution of transverse separations with the mean transverse 
size related to the inverse value of the mean transverse momentum of the quarks, 

{rl) = = 9 m P )) (33) 

where 



I\2 dkj* 


dW(t,z h ,k^,v)/d\n(k^) 


Ja 2 dkj< 


dW(t, Zh, kj,, v)jdk\ 





A qq dipole with starting value of transverse size r±(t p ) experiences evolution, 
forming the hadronic wave function and getting through the filtering process in the 
nuclear medium which absorbs dipoles of larger size more effectively. If it were in 
vacuum, it would take formation time tf which can be estimated as 

2z h u 

where h! is the first radial excitation of the hadron h. As was mentioned in the 
introduction, the formation time linearly rises with Zh, while the production time 
vanishes at z^ — > 1. 

The effective fragmentation function modified by the survival probability (trans- 
parency) of the pre-hadron in the nuclear matter reads, 

00 

D? ff {z h ,Q 2 ,v) = Jd 2 b J dzp A (b lZ ) 

—00 

x / dt W (t, z h ,Q 2 ,v)Tr(b,z + t, 00) . (36) 
Jo 

In the case of vanishing absorption, a^L — ► 0, this fragmentation function is identical 
to that of the free nucleon, except it is A times larger. The ratio of the effective 
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fragmentation function on a nuclear target to one on a free nucleon is identical to the 
experimental definition Eq. (5). 

As was mentioned, we use the dipole approximation, cr^(r) = Cr 2 . The factor 
C is proportional to the inverse mean hadron radius squared which is related to the 
parameters of the oscillatory potential, 

C=£f*- = « ro * w *«- (37) 

\ r T/h z 

where to = (m^ — mf l )/2 is the oscillatory frequency. For production of leading pions 
we took oo = 0.3 GeV and m q = 0.2 GeV which is constrained by the mean charge pion 
radius, 

(4) h = — = \ < r 2 > ch , (38) 
m q lo 3 

and nuclear shadowing data [43]. Using = 25 mb we arrive at C ~ 1.9. 

Note that the dipole cross section must be flavor independent, i.e. factor C must 
be the same for different species of mesons. This condition is satisfied by Eq. (37) 
as long as oc {rj)h- Such a relation between radiuses of hadrons and their cross 
sections is supported by data [44]. 



6 Comparison with data 

Unfortunately the low statistics of the HERMES data does not allow to perform a 
multi-dimensional comparison with the theory. One should integrate over all variables 
except one to get data with reasonable errors. 



6.1 Energy dependence of nuclear effects 

The results of numerical calculations for energy dependence of the ratio Ra{v) defined 
in (5) are depicted by dashed curve in Fig. 5 in comparison with with HERMES data. 
The cross sections for leading pion production were integrated over Zh and pt at the 
mean value of Q 2 same as in HERMES data. In what follows we always fix implicit 
variables at the values corresponding to the kinematics of HERMES 3 . 

Note that we ignored the induced energy loss in these calculations and integrated 
over zi- L > 0.5. Since the data are integrated over Zf- L > 0.2, this may cause some 
difference. We also neglected nuclear effects in data taken on a deuteron target when 
it is used as the denominator in (5). 

It is clear why nuclear suppression is maximal at low energies and has a tendency 
to vanish with increasing energy (data from EMC experiment [45] at much higher 
energies indeed demonstrate no nuclear effects). As we discussed in previous sections, 
hadronization process is completed by production of the leading colorless wave packet, 
which is then evolving into the final hadron over the formation time tf given by (35). 
If color neutralization happens outside the nucleus (t p > Ra), the evolution leads only 

3 We are grateful to Valeria Muccifora and Pasquale Di Nezza for providing us with required experimental 
information. 
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Figure 5: Comparison of the model predictions for the energy dependence 
of the nuclear ratio Eq. (5) for nitrogen and krypton targets HERMES 
data [2]. The solid and dashed curves correspond to calculations with 
and without the effects of induced radiation respectively 



to a phase shift, resulting in no nuclear effects, as long as we neglect the induced 
energy loss. 

On the other hand, the evolution may cause a strong attenuation if the wave packet 
is produced inside the nuclear target, as it happens at low energies where t p < Ra- 
At higher energies tp rises leaving less room for absorption. Besides, the formation 
time tf rises too, and the initial small size, ~ 1/Q, of the pre-hadron is evolving 
slower, leading to a lesser attenuation due to color transparency. Both effects lead to 
a stronger nuclear suppression at small than at high energies in accordance with data 
from HERMES [2] (Fig. 5), SLAC [46] and EMC [45,47] experiments. 

6.2 z^-dependence of nuclear effects 

Even at high energies the production time shrinks at large — ► 1 according to (8) and 
Fig. 2. Thus, at large Zh the path available for absorption of the pre-hadron reaches 
maximal length resulting in a strongest nuclear suppression. Such a shape of nuclear 
ratio falling toward Zh = 1 has been expected since 30 years ago [6,26], but high 
energy experiments could not reach sufficiently large values of Zh- Our predictions [1] 
are compared with EMC data in Fig. 6. We see that the interval of Zh where l p is 
short absorption is at work, is rather short and squeezed towards Zh = 1. This is 
because the energy of the EMC experiment is too high, but absorption is at work only 
if l p oc v(l—Zh) is shorter than the nuclear size, what needs quite a small (1— Zh) ocl/v. 

This is why measurements at much lower energies of HERMES were proposed 
in [1] in order to detect this effect. Fig. 7 demonstrates HERMES data for nitrogen 
and krypton in comparison with our predictions shown by dashed curves for the case 
of pure absorption. This is the first experimental confirmation for the longstanding 
prediction of vanishing production length at z^ — > 1. 

Note, however, that even at Zh — ► 1 nuclear suppression vanishes if both v and Q 2 
are sufficiently large. This is demanded in QCD by /^-factorization and is correctly 
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Figure 6: Comparison of the model calculations for the z^- dependence 
of nuclear suppression Rq u (5) at the energy < v >= 35 GeV (dashed 
curve) and < v >= 145 GeV (solid curve) with the EMC data [45,47]. 
Calculations are done at Q 2 = 11 GeV 2 . 



reproduced by our model. Indeed, at high Q 2 the initial size of the pre-hadron is very 
small, and it remains small long enough if the energy is high. Then color transparency 
eliminates nuclear effects. 



6.3 Q 2 dependence of nuclear effects 

At least two effects are sensitive to a variation of Q 2 , color transparency and the 
production length. Apparently, the initial size of the pre-hadron shrinks at larger Q 2 
and the nuclear matter becomes more transparent, i.e. the survival probability for the 
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Figure 7: Model predictions for the z h - dependence of nuclear ratios Ru N 
(5) and R»4 Kr calculated at values of < v > and Q 2 corresponding to the 
measurements at HERMES. The solid and dashed curve are the model 
predictions with and without inclusion of induced radiation, respectively. 



20 



pre-hadron increases. At the same time, the production length L contracts because of 
rising vacuum energy loss, and provides a longer path available for absorption. This 
effect leads to a stronger nuclear suppression at larger Q 2 . Thus, the two effects work 
in opposite directions and may partially or fully compensate each other. Indeed, it 
nearly happens: our predictions for the effect due purely to absorption depicted by 
dashed curves in Fig. 8 for neon demonstrate a rather weak Q 2 -dependence. Although 
the Q 2 -dependence demonstrated by the HERMES data is not steep either, one can 
see some difference. We will come to this point later after inclusion of corrections 
related to induced energy loss. 



7 Corrections for induced gluon radiation 

Multiple interactions of the quark during the time interval between the DIS point 
and production of the colorless pre-hadron cannot stop or absorb the quark, which 
keeps propagating through the medium, rotating in the color space. Nevertheless, the 
additional soft kicks gained by the quark force it to shake off more gluons. This addi- 
tional gluon radiation induced by the medium increases the loss of energy compared 
to vacuum. We calculate the induced energy loss in accordance with Eq. (3). This 
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Figure 8: Comparison of the model calculations for the Q 2 - dependence 
of the ratio Rn c of cross sections integrated over Zh with preliminary 
data from the HERMES experiment [3] . The solid and dashed curve are 
the model predictions with and without corrections for induced gluon 
radiation respectively. 
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result of [18] has corrections for finite medium and high energies [48,49]. However, the 
energies of HERMES are not high and the nuclear size is much longer than the gluon 
mean free path which is of the order of 0.2 — 0.3 fm [32,50,51]. 

Ai4 m Eq. (3) can be taken from data on p^-broadening in Drell-Yan reaction, or 
from calculations (see next section) which agree with data. Then one should add this 
induced energy loss to the vacuum one, i.e. to replace AE vac (t) =>■ AE vac {t)+AE in d{t) 
for the time-dependent quark energy, E q (t) = v — AE(t) in Eq. (21). 

Another modification of our calculations caused by p^-broadening in the medium 
is the modification of the hard scale. Namely, one should replace Q 2 =>- Q 2 + Ap?n(t) 
in Eq. (17) for the Sudakov factor and in Eq. (21). 

Both corrections enhance nuclear suppression, i.e. diminish the nuclear ratio. The 
results are shown by solid curves in Figs. 5, 7, 8 and 12. Apparently, these induced 
energy corrections vanish if l p — > 0. This happens at z^ — > 1 and at maximal Q 2 . 
Indeed, the dashed and solid curves in Figs. 7, 8 and 12 coincide in these limits. 

Note that inclusion of induced energy loss improves agreement with data. This fact 
may be considered as an indication to the importance of these corrections, although 
one should not probably expect very good description of data by such a rough model. 

8 Nuclear broadening of the transverse momen- 
tum distribution 

According to the perturbative QCD treatment of induced energy loss, it is caused 
by the broadening of transverse momentum of the quark initiating a jet during the 
hadronization process (at t < t p ). At longer time intervals, t > t p , the produced 
pre-hadron carries undisturbed information about the quark transverse momentum 
(since it has no interactions) which thus can be directly measured. Comparison with 
theoretical predictions seems to be of special importance. 

A parton propagating through a medium experiences multiple interactions in the 
medium and performs Brownian motion in transverse momentum plane. As a re- 
sult, the mean transverse momentum squared of the quark linearly rises with the 
path length. Quantitatively this process has been successfully described within the 
light-cone dipole approach [19, 20] (otherwise, most of models just fit the data to be 
explained). If the DIS takes place on a nucleon with coordinates (b, z), the medium 
modified transverse momentum distribution is expressed in terms of the universal cross 
section of a qq dipole with a nucleon, 



dN^(b, z) 



q 



(PpT 

x exp 



J d 2 rx d 2 r' T Vi q {rT, r* T ) exp [ipxi^r — t't) 



T A (b,z,z + l p ) N _, 

5 VqqilT ~ r T , V) 



(39) 



Here T A {b, z,z + l p ) = j* +lp dz' PA(b,z') is the nuclear thickness function; cr^(rT, f) 
is the phenomenological dipole cross section fitted to data for the proton structure 
function and photoabsorption cross section. It includes nonperturbative effects, as well 
as all corrections for gluon radiation (which give rise to the energy dependence of the 
cross section) also contributing to the broadening of the quark transverse momentum. 
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The quark density matrix in coordinate representation, 



k 2 

Q, q (r T ,r' T ) = — exp 

7T 



f(4 + ^ 



(40) 



where fco is the mean transverse momentum of the quarks, describes the initial dis- 
tribution of the quark in the target nucleon. It controls the transverse momentum 
distribution of the valence quark emerging from DIS on a free nucleon. Therefore, the 
mean transverse momentum is of the order of the inverse proton size, i.e. is rather 
small, ko ~ ^-qcd- Usually data on hard processes need a much larger primordial 
transverse momentum. This effect is related to the next-to-leading order corrections 
related to hard gluon emission. Since we explicitly include gluon radiation in our 
model, additional next-to-leading order corrections would lead to double counting. 

Note that our parameter-free approach successfully explained available data for 
the Cronin effect in pA collisions in fixed target experiments, and correctly predicted 
it for d — Au collisions at RHIC [52, 53]. Data for nuclear broadening in Drell-Yan 
reaction also are well reproduced [20]. 

To get the hadron transverse momentum distribution we should convolute the 
quark p^-distribution, Eq. (39) calculated without the exponential factor, with the 
fragmentation function D h / q (zh,PT), 



a DIS d Pt J d IT 

dD h/q (z h ,k 2 ) ( k T \ 
x ^ 6 ^pr — (IT —J . (41) 

Now we can make use of our model for the vacuum fragmentation function D h / q (zh, Q 2 ), 
Eq. (25), and calculate the mean momentum squared of the produced meson. It hardly 
varies with Q 2 within the kinematic range corresponding to HERMES data. We found 
(p*) = {k^)/A + kl = 0.3 GeV 2 , where k = 0.2 GeV. This result agrees with the 
pr-distribution of hadrons measured at SLAC [54] at energies overlapping with the 
HERMES energy range. 

The pT-dependence of the effective fragmentation function in the nuclear medium, 
according to Eq. (36), reads, 



7 UTS d PT J J 



'DIS 

— oc 



r°° , s f o o dN A (b,z) 

/ dt Tr(b, z + t,oo) / d 2 k T d 2 q T — V 
Jo J d A q 



dW(t,z h ,k^,u) ^ _ k T \ 
x ^- S[p T -q T - Y j (42) 

We also included the small effect of Fermi motion in the nucleus. 

The calculated ratio of the two distributions, nucleus-to-nucleon, is depicted in 
Fig. 9 for nitrogen and krypton in comparison with data from the HERMES experiment 
[2]. The model reproduces well the data with no adjustment, although the curves 
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Figure 9: The pr-dependent nucleus-to- nucleoli ratios of the cross sec- 
tions of hadron production in DIS on nitrogen (solid curve) and krypton 
(dashed curve) in comparison with HERMES data. 



seem to rise somewhat steeper than the data at large pr- That may be caused by the 
difference in the bottom limits of integration over (0.2 in the data, versus 0.5 in 
our calculations). 

9 Disentangling absorption and induced energy 
loss 

The production length is a fundamental characteristics of the early stage of hadroniza- 
tion, and it is important to sort out those observables which are sensitive to l p . 

Such observables could also help to disentangle different models of hadronization. 
The current estimates of l p within different models (strings [2,6,10,26], pQCD [1]) 
are pretty close and show that l p is rather short compared to the nuclear sizes for the 
kinematics of the experiments at HERMES [2] and especially at Jefferson Lab [55]. 
At the same time, some models [4] make an assumption that color neutralization 
happens always outside the nucleus. Although this is an ad hoc assumption with no 
justification, would be useful to confront it with direct experimental evidence. 

9.1 Flavor dependence 

The survival probability for the produced pre-hadron depends on the wave function of 
the final hadron, as one can see in Eq. (27). Also intuition tells us that the projection 
to a wave function with a larger mean radius should emphasizes larger qq separations 
in the pre-hadron, leading to a stronger absorption. Since the pion has a larger 
radius than the kaon, the former is expected to be produced with a stronger nuclear 
suppression than the latter. This may be a good signature for absorptive effects and 
shortness of l p . On the contrary, in the case of l p 3> Ra there should be no sensitivity 
to the hadronic size. The two curves are different only by the value of the charge 
radius, {r 2 ) c h (kaon or pion) used in calculations. 
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Figure 10: Comparison of nuclear effects for kaon (thin curve) and pion 
(thick curve) productions as function of energy in comparison with HER- 
MES data. 

The HERMES experiment has recently provided high quality data for hadron pro- 
duction with particle identification [2]. Here we concentrate only on production of 
pions and positive kaons, since negative kaons and antiprotons do not contain valence 
quarks common to the target nucleon, while proton production may be contaminated 
with the target nucleons. 

Our predictions for the production rates of positive kaons and pions integrated 
over Zh are depicted in Fig. 10 as function of photon energy together with recent data 
from HERMES. 

Data confirm the expected identical suppression for positive and negative pions, 
and less attenuation for kaons. Although our calculations underestimate kaons, this 
is due to the difference in the bottom limit of integration over z^, {zhjmin = 0.2 in 
the data, and 0.5 in our calculations. At small Zh kaons are copiously produced via 
reaction irp — ► KA which has a very low energy threshold. It would make sense if 
experimental data used a larger value of (zh) m in- 

Although the alternative model, the energy loss scenario, is not sensitive to the 
hadronic radius, one may expect a difference in quark fragmentation functions for 
q — > 7T and q — > K, where q is a light quark, u or d. Indeed, the end-point behavior of 
the fragmentation function, at Zh — > 1, is dictated by the Regge phenomenology [31], 

D v/q (z h ) oc (l-* ft )i-«p(o) = (i_^)0.5. 
D K+/q (z h ) oc (1-^) 1 -^(0) = (1-^) - 8 . (43) 

Since the induced energy loss leads to a shift in Zh => Zh + Az/,, it should lead to 
a stronger suppression of kaons which have a steeper fragmentation function. This 
expectation is at variance with HERMES data depicted in Fig. 11. 

At the same time, the data agree well with our predictions shown by thick and thin 
curves for pions and positive kaons respectively. This fact confirms our claim that the 
disagreement with data for kaons seen in Fig. 10 is due to presence of small- Zh events 
in the data. 
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Figure 11: Nuclear ratios for pions (thick curve) and positive kaons (thin 
curve) as function of z h predicted by our model in comparison with HER- 
MES data. 



Note that the maximum effect of flavor dependence should be expected for shortest 
l p , e.g. at Zh — > 1. Indeed, our calculations predict quite a large difference between 
nuclear effects for K + and pions at Zh — > 1, as is demonstrated in Fig. 11. Unfortu- 
nately, the accuracy of the data is not sufficient to see such a variation of the flavor 
dependence with Zh- 

9.2 Q 2 -dependence 

The fit to HERMES data performed in Ref. [4] within the energy loss scenario led 
to the value of the universal parameter, twist-4 quark-gluon correlation tensor, which 
is twice as small as extracted from Drell-Yan data and 25 times smaller than follows 
from data on high-py dijet production. Such a steep variation of the parameter which 
must be process-independent was attributed in [4] to a strong scale dependence. If it 
were true, the induced energy loss effect would rise with Q 2 . Accordingly, the nuclear 
suppression should become stronger at larger Q 2 . This expectation is in obvious 
contradiction with the HERMES results depicted in Fig. 8. 
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Figure 12: The same as in Fig. 8, but predicted for krypton. 
On the other hand, one may expect a more pronounced rise of the nuclear ratio 
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with Q 2 for heavy nuclei. Our results for krypton depicted in Fig. 12 confirm this 
expectation, although the difference from data on neon target, Fig. 8 is not a dramatic 
effect. Note the sizeable variation of nuclear suppression with Q 2 at low energies. In 
this case the production time is short at any Q 2 and the variation of nuclear effects 
versus Q 2 is only due to color transparency. This is why the nuclear ratio rises with 
Q 2 even without induced energy loss corrections. 

9.3 Zh dependence of p^-broadening 

According to the discussion in the previous section, nuclear broadening of the trans- 
verse momentum distribution might be the most sensitive probe for the production 
length, and provides a direct measurement of l p , since Ap 2 ^ oc l p . Indeed, only the 
hadronizing quark at short time intervals t < t p contributes to the broadening. As 
soon as the pre-hadron is created, no further broadening occurs, since inelastic inter- 
actions are prohibited for the pre-hadron, and only broadening via elastic rescattering 
is possible. However, the elastic cross section is so small that even for pions the mean 
free path in nuclear matter is about 20 fm. It is even longer for a small-size pre-hadron 
due to color transparency. Therefore, we should expect a disappearance of the broad- 
ening effect at large Zh —>■ 1 since t p — ► 0. In Fig. 13 we show our predictions for the 
nuclear ratio as function of pr for different Zh bins. 
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Figure 13: p^-dependence of nucleus-to-nucleon ratios binned in Zh for 
hadroproduction in DIS on nitrogen (left) and krypton (right). 

Note, however, that broadening should be weaker at larger Zh for any model which 
is able to describe the observed z/i-dependence of the nuclear ratio (see Fig. 7). Indeed, 
the fact that RA(zh) < 1 implies that DIS happens only in a part of the nuclear volume, 
namely its back side. Then, the fact that R^z^) decreases towards = 1 means that 
a smaller part of the nuclear volume is working at larger Zh, i.e. the DIS reaction is 
pushed towards the back surface of the nucleus. Therefore, even in the energy loss 
model the path available for broadening (from the DIS point and on) becomes shorter 
at large zy t . This should also lead to a reduction of the broadening, but rather small, 
only about 10% for nitrogen and 20% for krypton. At the same time, in our approach 
the broadening completely vanishes at z^ — > 1 since l p — ► 0. It would be extremely 
important for HERMES and the experiment running at Jefferson Lab [55] to provide 
relevant data. 
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9.4 Double hadron production 



One may think that a process of double hadron production is a sensitive tool which is 
able to disentangle the inside-outside hadronization. If absorption of the pre-hadrons is 
the main contributor to nuclear suppression, naively one would expect the nuclear ratio 
for double-hadron production to be squared, i.e. quite small, compared to inclusive 
single-hadron channel. 

A closer look, however, shows that the situation is quite complicated and such 
an expectation is oversimplified. Let us assume that the two pre-hadrons are created 
inside the nuclear medium with production lengths l\ and I2 respectively, and I2 > h- 
For the sake of simplicity, just for this consideration, we assume that the nuclear 
density is constant, both hadrons have the same constant absorption cross sections a, 
and we neglect any additional suppression related to induced energy loss. Then it is 
straightforward to write the nuclear suppression ratio at given impact parameter b, 

00 00 
R^ h \b) = J dz pa(P, z) exp —a J dz'pA(b,z') 

-oo z+h 
oo 

x exp -a J dz' p A (b,z') (44) 

Using the approximation of a constant nuclear density, pA^r) = paQ(Ra — r), and 
neglecting exp(— 2apAL) <C 1, where L = ^ R 2 A — b 2 , we arrive at the relative, double 
to single hadron production, nuclear suppressions, 

R ( A h \ b ) eM-ap A Al) + up A h 



R^ib)' l + <?PAl P 



(45) 



where Al = I2 — l\ and l p is the single-hadron production length. This ratio varies 
dependent on the values of l\, I2, l p and the mean free path, 1/apA, in the medium. 
In the case of a very dense matter, e.g. produced in heavy ion collisions, one should 
expect lop a 3> 1. Contrary to naive expectations, the production rates of single and 
double hadrons according to (45) should be suppressed equally. The data for high- 
PT production at RHIC (see next section) seem to confirm this. This situation may 
change at very large px due to contraction of the production lengths [see Eq.(48)]. 

In the case of DIS in the dilute nuclear medium production lengths are comparable 
with the mean free path, which is about 2 fm. Then the result should vary dependent 
on values of Zh for each of the hadrons. For instance, if values of Zh in both cases 
approach their maximal values, i.e. z^ — > 1 and Zh\ + z^ — > 1, all production lengths 
vanish, and the ratio Eq. (45) equals 1/2, i.e. a double-hadron is suppressed twice as 
much as a single-hadron. 



10 Heavy ion collisions: Energy loss or absorp- 
tion? 

The situation in the process of high-p^ hadron production in heavy ion collisions at 
very high energies is similar to the process of hadron production in DIS off nuclei at 
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medium energies. Indeed the cartoon in Fig. 14 shows the development of hadroniza- 
tion after the nuclear disks passed through each other. Partons are copiously created 
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Figure 14: Propagation and hadronization of a high-fc^ parton in the 
dense partonic matter created right after collision of two nuclei in their 
cm. frame. 



from vacuum forming a dense medium of about the transverse size of the nuclei. The 
high-pr parton which was created in a hard partonic reaction at the earliest stage of 
nuclear collision is propagating through this medium. This process ends up with pro- 
duction of a hadron detected at macroscopic distances. The momenta of the hadrons 
in the cm. frame of collisions ranges from a few up to 10 GeV (in the available 
data), which overlaps with the kinematics of HERMES [2] and experiment running at 
Jefferson Lab [55]. 

The kinematics and interpretation of this process are much less certain than in 
DIS. In the latter case the initial quark energy is known, as well as z^, and the density 
and size of the medium. On the contrary, the production cross section of a high-fcj- 
parton is a result of convolution of the initial parton distributions, hard partonic cross 
section, and fragmentation functions. Therefore the initial parton energy and values 
Zh are not known. Besides, one has to sum over different species of partons. 

These measurements are aimed at probing the density of the created medium, 
which is varying in space and time. This brings even more model dependence to the 
interpretation of data. In view of these uncertainties in interpretation of the data on 
heavy ion collisions, one desperately needs a reliable model of in-medium hadronization 
tested in more certain situations like DIS. 

There is a principal difference between hadronization processes in DIS and high-pj- 
hadron production at mid rapidities. The dependence of the production length, l p , 
on the jet energy varies substantially depending on the process. Indeed, our previous 
experience says that the production time is proportional to the jet energy (the Lorentz 
factor). In soft interactions, 

La-. (46) 

K 

The string tension k in the denominator provides the correct dimension for L. 

In the case of DIS the production time is also proportional to energy, but depends 
on the photon virtuality, 

l P oc . (47) 
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Here the right dimension is provided by the photon virtuality, Q 2 . The physics is 
transparent: the higher is Q 2 , the stronger is the kick from the photon, the more 
intensive is gluon radiation and vacuum energy loss by the quark. Correspondingly, 
the shorter must be the color neutralization time, l p , if one wants to produce a leading 
hadron with large Zh- 

Now we are coming to the point. In the case of 90° parton scattering, like in 
the process illustrated in Fig. 14, the jet energy [z/(DIS)=>- &t(RHI)] and the parton 
virtuality [Q 2 (DIS)=^ fc^(RHI)] are controlled by the same parameter, which is the 
transverse momentum of the parton, kx- Therefore, in this case, 

l p oc -i- . (48) 

It is clear why the energy dependence of the production time switches to the inverse. 
The vacuum energy loss with a rate proportional to k T is so intensive, that in spite of 
the Lorentz factor kx, the hadronization process must finish shortly after the hard par- 
tonic collision, otherwise the parton energy will degrade too much, making impossible 
production of an energetic hadron. 

It is also clear that in pure perturbative QCD calculations the dimensional param- 
eter, Aqcd, can emerge only under a log. Therefore, Eq. (48) presents the unique 
possibility to provide the correct dimension for l p . Thus, the higher is the energy of 
the jet, the faster the leading pre-hadron is produced. This outrageous conclusion may 
contradict simple intuition based on the experience with the string model. 

One may think that no pre-hadron can be produced inside a deconfined medium. 
This is not correct. Since the mean radius of the produced pre-hadron decreases like 
1/kx, it should be smaller that the Debye screening radius at large kx- Later on, the 
pre-hadron is evolving its size and may be dissolved in the medium, but that would 
just mean a large absorption cross section. 

The energy dependence of the production length suggested by Eq. (48) for high-pr 
processes is inverse to what one has for DIS, Eq. (47). Correspondingly, the energy 
dependence of nuclear suppression should be opposite to DIS. As long as the nuclear 
ratio in DIS rises with energy (Figs. 5-10), we should expect a falling p^-dependence 
of the nuclear ratio for heavy ion collisions. This unusual effect was indeed observed 
at RHIC. Data from the PHENIX experiment depicted in Fig. 15 demonstrate quite 
a strong fall of the ratio at px > 2 GeV. Of course at high px one should expect 
saturation (l p = 0) and the following rise due to color transparency. 

On the contrary, in the energy loss scenario, there is no principal difference between 
the two cases. As long as the nuclear ratio rises with energy in DIS, it must rise with 
Pt in heavy ion collisions (unless one makes different models of hadronization in the 
two reactions). Such a behavior was indeed predicted in [4] as is shown in Fig. 15. 
That was a true prediction before the data at high px were released. Of course, once 
the data are known, one can find an explanation by introducing new, exotic physical 
assumptions. 

11 Conclusions and discussion 

The model of in-medium hadronization [1] developed long ago (last century) has 
successfully predicted the nuclear effects for inclusive electroproduction of leading 
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Figure 15: Data from the PHENIX experiment [56] for pion suppression 
in 10% central gold-gold collisions at y/s = 130 GeV in comparison with 
the prediction from [4]. The figure is borrowed from Ref. [56]. 



hadrons. It was realized that the string model is rather irrelevant to the early stage of 
hadronization in such a perturbative QCD process as DIS. Instead, there has been de- 
veloped a perturbative description which possessed a clear pattern for the space-time 
evolution, and allowed to carry out numerical estimates. The key points of the model 
are as follows. 

• A leading quark originated from DIS loses energy for hadronization which we 
treat perturbatively via gluon bremsstrahlung. It is important to discriminate 
between vacuum and induced energy losses (which are frequently mixed up). The 
former is always present, even on a free nucleon target, and is the main source 
of energy loss. In the case of in-medium hadronization the quark radiates more 
gluons due to multiple soft interaction with the medium. This correction is 
usually much smaller than the vacuum energy loss. 

• It is important to discriminate between production and formation times (which 
are frequently mixed up). Energy loss (both vacuum and induced) stops when 
color neutralization occurs, i.e. the leading quark picks up an antiquark and 
produces a pre-hadron, which is a colorless qq dipole with no stationary wave 
function. The corresponding time interval is called the production time, or pro- 
duction length. Energy conservation enforces the production time to vanish 
at the kinematic limit of production of hadrons with maximal possible energy 
(z h - 1). 

• The pre-hadron attenuates on its way out of the nucleus with an absorptive 
cross section which is controlled by the varying dipole size. Color transparency 
is an important ingredient of this dynamics. Eventually, on a longer time scale, 
called formation time (or length) the dipole develops the hadronic wave function. 
We employ a rigorous quantum-mechanical description of this process within the 
light-cone Green function approach. Contrary to the production length vanishing 
at Zh — ► 1, the formation length reaches a maximal value in this limit. 

• All these effects important for nuclear modification of the hadron production 
rate can be evaluated. Then the model is able to predict in a parameter free 
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way the nuclear effects in DIS as function of energy, Zh, Q 2 and px- Data from 
HERMES nicely confirm the predictions. 

• Searching for observables which are able to disentangle different models, we con- 
cluded that the energy loss scenario should have problems explaining the already 
available data for flavor and Q 2 dependence of nuclear suppression. This model 
is insensitive to the hadronic size, but the difference in fragmentation function 
at Zh — > 1 leads to more suppression for kaons than for pions. This is at vari- 
ance with HERMES data. Also the comparison of data on different reactions at 
different hard scales led to a conclusion that pr broadening and energy loss rise 
with Q 2 . Such an expectation also contradicts HERMES data. It would be very 
informative to have data on the variation of the p^-distribution with z^ , which 
is a direct way to measure the production length. 

• The new feature of hadronization dynamics in high-p^ processes is an inverse 
dependence of the production time on the jet energy compared to what is known 
for DIS. Correspondingly, the pr-dependence of nuclear suppression in heavy ion 
collisions should be enhanced at larger px, contrary to DIS where nuclear effects 
vanish at high energy. Such a puzzling behavior was indeed disclosed by recent 
measurements at RHIC. Nevertheless, at higher px we expect the effect of color 
transparency to take over, then the £>r-factorization will be restored. 
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